Isentropic Curves at Magnetic Phase Transitions 
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Experiments on cold atom systems in which a lattice potential is ramped up on a confined cloud 
have raised intriguing questions about how the temperature varies along isentropic curves, and how 
these curves intersect features in the phase diagram. In this paper, we study the isentropic curves of 
two models of magnetic phase transitions- the classical Blume-Capel Model (BCM) and the Fermi 
Hubbard Model (FHM). Both Mean Field Theory (MFT) and Monte Carlo (MC) methods are used. 
The isentropic curves of the BCM generally run parallel to the phase boundary in the Ising regime of 
low vacancy density, but intersect the phase boundary when the magnetic transition is mainly driven 
by a proliferation of vacancies. Adiabatic heating occurs in moving away from the phase boundary. 
The isentropes of the half-filled FHM have a relatively simple structure, running parallel to the 
temperature axis in the paramagnetic phase, and then curving upwards as the antiferromagnetic 
transition occurs. However, in the doped case, where two magnetic phase boundaries are crossed, 
the isentrope topology is considerably more complex. 

PACS numbers: 05.10.Ln, 71.10.Fd, 75.30.Kz,75.10.Jm,71.30.+h,64.60.Cn 



I. Introduction 

In systems of interacting degrees of freedom, decreas- 
ing the thermal fluctuations often leads to the formation 
of ordered states. The traditional, and natural, mea- 
sure of these fluctuations is the temperature T itself, 
which then forms one axis of the associated phase dia- 
grams. However entropy S can also be used to quantify 
the amount of disorder. Indeed, a phase diagram using 
S as an axis naturally provides a somewhat different per- 
spective on the topology of the ordered and disordered 
regions- since the entropy changes more rapidly where 
transitions occur, it magnifies these interesting portions 
of the phase diagram. 

Recent experiments^— on trapped ultracold atoms in 
optical lattices have provided a further motivation for 
employing the entropy as one of the variables in describ- 
ing phase diagramsj££ In these systems, the temperature 
and entropy of the atomic cloud are known prior to the 
adiabatic ramp-up of the optical lattice, but the precise 
change in temperature during this process is uncertain. 
Thus the determination of the entropy values at which 
various phenomena occur, like local moment formation, 
magnetic ordering, and so forth, is important, supple- 
menting the more typical discussion of the temperatures 
at which these phenomena occur. Furthermore, since at- 
taining low temperatures is crucial for the emulation of 
many-body ordering effects seen in solid state systems, 
a central question is whether T rises or falls (adiabatic 
heating or cooling) as the optical lattice is established. 

This work examines the isentropic curves of two mod- 
els of magnetic phase transitions. The two-dimensional 
Blume-Capel Model (BCM)I£ is studied first, and its 
thermodynamics are computed both in Mean Field The- 
ory (MFT) and with Monte Carlo (MC) methods. An 
itinerant (quantum) Hamiltonian, the Fermion Hubbard 
Model (FHM) is next examined within MFT. The isen- 
tropes of the FHM are compared with recent Quantum 
Monte Carlo (QMC)2 results at half-filling, a density at 



which QMC can be performed to low temperatures. The 
isentropes are also computed when the system is doped 
away from half-filling, a parameter regime where phases 
with long range order are inaccessible to QMC. The BCM 
and FHM form an interesting pair of models to compare, 
since both contain two energy scales, one which controls 
the density and the other which tunes the strength of 
spin-spin interactions. 

The paper is organized as follows: Sec. II presents both 
the BCM and the FHM, and the computational methods 
used. Sec. Ill then details the isentropes of the BCM, 
while Sec. IV focused on the FHM. Sec. V summarizes 
and further discusses the results. 



II. Models and Computational Approach 



A. Blume Capel Model 

The Blume-Capel Mode££ is, 

E = -j^2SiSj+ D^S? . 

<ij> i 



(1) 



The spins Si can take three values Si = 0, ±1, where 
the value Si = can be regarded as a 'vacancy'. The 
first term in the energy represents a ferromagnetic cou- 
pling J between near-neighbor spins. The second term 
provides an 'impurity' chemical potential D for the va- 
cancies. When D is large and negative, vacancies are 
suppressed, and the BCM maps onto the Ising model. 
Here we will consider a square lattice geometry. 

The presence of a three component spin gives rise to 
the possibility of first order transitions and tricritical be- 
havior, as first emphasized in [Toj . Part of the original 
motivation for the BCM was to provide a description of 
the tricritical phenomena induced by 3 He vacancies in 
supcrfluid 4 HeAi~— Since its introduction, the BCM has 
been extensively studied^ - — both in the form given in 
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Eq. [T] and in several variants which include additional 
terms in the energ y 22 ' 23 and generalizations to vector 
spins which more correctly capture the continuous sym- 
metry of the supcrfluid order parameter j2£^i 

The solution of the BCM within MFT is straightfor- 
ward. The results for the free energy and entropy (per 
site) are, 

/ = - jj In [2 cosh(/3 Jzm) + e 0D ] + D (2) 

s = -l3[Jzm 2 + (1 - p)D] + ln[2 cosh(/3 Jzm) + e pD } 

where z is the coordination number (z = 4 for square 
lattice). The magnetization m and density of spin states 
p are given by the self-consistency equations: 

2 sinh(/3 Jmz) 

TTL == 

e? D + 2 cosh(/3 Jmz) 
2 cosh(/3 Jmz) 
9 = ef D + 2 cosh(/3 Jmz) ( ' 

Here we set J = 1 as a unit of energy, and solve Eqs. 12131 
to obtain the free energy and entropy for any values of 
independent variables D and . 

We compared two different MC approaches for com- 
puting entropy in the BCM: thermodynamic integration 
and the Wang-Landau (WL) algorithm. With thermo- 
dynamic integration, the entropy is computed by making 
multiple MC runs at different, fixed inverse temperatures 
/? and integrating, 

s(/3) = S (oo)+/3e- / e((3)d(3 , (4) 

where s and e are the entropy and energy per site, re- 
spectively. Equation |4] is obtained by applying integra- 
tion by part a 9 ' 28 ' 29 to the standard relation of entropy 
to the specific heat: S(T) = J Q T C{T')/T' dT'. For the 
BCM, s(oo) = In 3, reflecting the three possible choices 
of spin. In most instances, this approach produces bet- 
ter results, since it does not rely on the determination of 
specific heat C(T), which is noisier than the energy e(T). 

The WL algorithm is a flat-histogram MC method for 
calculating the density of states g(E). In this method, 
a random walk is performed in the energy space of the 
BCM, sampling E with a probability proportional to 
l/g(E) and adjusting the distribution of g(E) until each 
energy (E) value has close to the same probability. Ul- 
timately, this process produces a flat histogram of oc- 
currence for all energy states in the random walk. Since 
the density of energy states that results is independent 
of temperature, we can compute the partition function 
Z = g(E)e~P E for any temperature. Consequently, 
we can determine the value for any thermodynamic vari- 
able of interest- in our case, the free- energy and entropy- 
without performing multiple MC simulations at different 
0- 

Figures and analysis reported for BCM will be those 
from WL results on 16x16 lattices. The values obtained 



for the entropy with the two methods, however, were 
found to be equivalent to within a fraction of a percent. 

B. Hubbard Hamiltonian 

The Fcrmion Hubbard Model ^ 

OS)* 

i icr 

describes the magnetism of itinerant electrons, in con- 
trast to the static spins of the Blume-Capel model. Here 
ci _(c. ) are creation(destruction) operators for fermions 
of spin a on lattice site i, and n\ a are the associated 
number operators, t = 1 sets the kinetic energy scale 
for the hopping of fermions between near neighbor sites 
of a square lattice. U is an on-site energy cost for dou- 
ble occupancy, and the chemical potential (i controls the 
filling. 

One interesting property of the square lattice, near- 
neighbor hopping is that the associated dispersion rela- 
tion e(k Xl k y ) = — 2t[cosk x + cosk y ] has a (logarithmic) 
divergent density of states at half-filling. As a conse- 
quence, the Stoner criterion suggests that an arbitrarily 
small interaction U will induce a magnetic instability. 
This is reflected in the phase diagrams shown in Figs. [5] 
and [5] 

The FHM has been widely used to study strong cor- 
relation effects in solids, from magnetism to metal- 
(Mott) insulator transitions, and high temperature 
supcrconductivity^2r— Recently, the FHM and its 
bosonic counterpart have attracted considerable inter- 
est for describing the behavior of cold atoms trapped in 
an optical lattice produced by interfering laser beamsi 
Compared with traditional condensed matter experi- 
ments, these optical lattice systems are thought to more 
precisely mimic the FHM, while allowing tunable con- 
trol over parameters like U and t. As mentioned in the 
introduction, this provides a central motivation for this 
paper. 

Our MFT approach is the usual one in which each term 
of the interaction is decoupled: Un^n-^ — > Un^fa^) + 
U(ntf)ri\± — U(ni^)(riii). The resulting quadratic Hamil- 
tonians H a arc diagonalized, and the expectation val- 
ues (riia) recomputed. The process is iterated to self- 
consistency, a process which minimizes the free energy 
F. To be somewhat more precise, the MFT calculation 
is actually performed in momentum space by making par- 
ticular, simple, paramagnetic (P), ferromagnetic (F), and 
antiferromagnetic (AF) ansatz for the real space expecta- 
tion values (riia). Such a choice does not allow for more 
complex spatial patterns of charge and spin such as are 
present in striped phases^ 

Our results for the isentropic curves of the FHM will be 
obtained exclusively within MFT. It is possible to obtain 
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these curves using more exact approaches like QMC^ but 
only in parameter regimes like half-filling where there is 
no sign problem. Previous Dynamical Mean Field Theory 
work has also reported data for the isentropes on a cubic 
latticed 



III. Isentropic Curves of the Blume-Capel Model 

Single Site Limit 

We begin our analysis of the isentropic curves by dis- 
cussing the J = limit of a collection of independent 
spin-1 sites. The partition function is, 

Z=l + 2e~ pD , (6) 

from which we can derive the internal energy, free energy, 
and entropy, 

(E) = Z'^De'^ 

f = -imz 

S = f3((E)-F). (7) 

Since only the combination /3D enters the expressions 
in Eq. [Jj the isentropic condition, S = So implies that 
(3D = c, where the constant c is obtained by solving the 
transcendental equation, 

So = -^^ + ln(l + 2 e - c ). ^ 

1 + Ze c 

We immediately see that the isentropic curves are 
straight lines D = cT. Note that for values of So greater 
than In 2 ~ .693, there are two solutions for c = D/T, 
one positive and one negative, representing positive and 
negative values of D. In this single site limit, values of 
D less than zero will suppress all vacancies, preventing 
entropy values less than In 2. 

Despite the simple nature of this calculation, it allows 
us to make some immediate statements about adiabatic 
heating and cooling. We see that in traveling along an 
isentrope of increasing D from D = — oo, the temperature 
T decreases, since one moves along a line with c < 0. 
Ultimately, one reaches D = where the sign of c changes 
to positive, and further movement along the isentrope 
results in an increase in temperature. 

The interesting question is, of course, how the topology 
of the isentropes of the noninteracting systems is altered 
by interactions and, in particular by the strong collective 
effects which occur near the magnetic phase boundary. 

Mean Field Theory 

Having discussed the J = limit, we now address the 
isentropic curves of the BCM treated within MFT. Solv- 
ing the self-consistency Eq|31 we show the resulting isen- 
tropes in Fig.Q]along with the mean- field phase boundary 
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FIG. 1: (color online). Mean-Field Theory calculation of the 
isentropes of the Blume-Capel model. The temperature T 
decreases as D increases along the isentropes in the ordered 
phase, and then rises in the paramagnetic phase. S jumps at 
the phase boundary below the tricritical point, and is contin- 
uous above. 

(PB). The PB consists of a first order line (shown in full 
red) and second order boundary (dashed green), which 
meet at the tricritical point (TP): (T tp , D tp ) = (§, § In 2). 
The TP and second order PB can be obtained analyti- 
cally by expanding the free energy in powers of the mag- 
netization m and computing the Landau coefficient for 
the to 2 term, thus fixing the transition line. (See Ref. 0.) 

At the first order phase boundary, the entropy shows a 
characteristic "jump," or discontinuity, as the isentropes 
take a sharp jog to the left (smaller T) in passing through 
the boundary (Fig. [I}. The entropy decreases in this 
traversal of the PB to higher D values at fixed T due to 
the large reduction in the spin density p. Below the PB 
in the ferromagnetic phase (F), spin-0 "vacancies" are 
sparse, resulting in an effective (Ising-like) two spin state 
region and lower entropy. Above the PB, all three spin 
states are present with higher entropy. 

In contrast, at the second order boundary (dashed 
green line in Fig. [1}, the isentropes are continuous, but 
exhibit a change in slope at the PB. As we will see shortly, 
the precise details of this behavior are peculiar to the 
MFT and do not carry over to the exact (MC) solutions 
which incorporate short range correlations and fluctua- 
tions. Nevertheless the general topology of the isentropes 
in MFT agree with those of MC. 



Monte Carlo- Wang-Landau 

The phase diagram of the BCM was determined by 
first using the WL density of states to compute thermo- 
dynamic variables, the free energy /, entropy s, specific 
heat C(T), and magnetic susceptibility x(T), for a fine 
grid of points in the (£>,T)-planc. We then examined 
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FIG. 2: (color online). Wang-Landau MC results for the spe- 
cific Heat C(T) for different values of the vacancy chemical 
potential D = 0.0, 0.5, 1.0, 1.5 in the second order phase tran- 
sition region. The peak positions give one estimate of the 
phase boundary. 



FIG. 3: (color online). Wang-Landau MC results for the spin 
density contours of the Blume-Capel model. The vacancy 
density increases as the temperature T or chemical potential 
D rise. Ferromagnetic order is lost when roughly one third of 
the sites are vacant. 



specific conditions to locate the phase transition points. 
For 1st order transitions, we looked for jumps in the en- 
tropy coinciding with a vanishing magnetic order param- 
eter (m). For 2nd order transitions, we located peaks in 
specific heat and magnetic susceptibility curves plotted 
as a function of temperature for fixed D. 

Representative results arc shown in Fig. [2J Here C(T) 
is obtained for fixed impurity chemical potentials D = 
0.0, 0.5, 1.0, and 1.5. Peaks are observed at temperatures 
which are associated with the large fluctuations at the 
second order phase transition. Finite size scaling, and 
other sophisticated techniques can be exploited to locate 
T c very precisely. For example, as D is increased further, 
it is known that a tricritical point exists in the BCM at 
(T t p,£>tp) = (.610,1.965). This change in behavior can 
be seen numerically by monitoring, the appearance of 
hysteresis loops when MC simulations arc done sweeping 
D at fixed T <T tp . 

Plots like that of Fig. [5] for different values of D, as 
well as sweeps in D at fixed T, were used to locate the 
phase transition line in the T — D plane. This is shown 
as the black dashed line in Fig. 2J The phase boundary 
thus obtained is in good agreement (less than a percent 
difference) with published results. The new aspect of 
Fig. Uis the inclusion of the isentropic curves. 

The gross linear structure of the isentropes for large 
T I J and D / J is well explained by the J = (single site) 
analysis earlier in this section. However, at intermediate 
T/J and D/J, the isentropes are rounded, especially in 
the vicinity of the phase boundary. Indeed, the boundary 
between adiabatic heating and cooling docs not precisely 
follow the PB but instead occurs along a separate trajec- 
tory of somewhat larger impurity chemical potential. 

As discussed in (37}, the key qualitative feature of the 
isentropes is that they move out to higher T as they leave 
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FIG. 4: (color online). Wang-Landau MC results for the isen- 
tropes of the Blume-Capel model. As for the MFT calcula- 
tion, the isentropes exhibit adiabatic cooling in the ferromag- 
netic phase, and heating in the paramagnetic phase. The 
entropy curves become discontinuous at the tricritical point. 

either side of the phase boundary (into the paramagnetic 
or ferromagnetic phases). The simple picture of this re- 
sult is that the boundary represents a line of a high de- 
gree of competition between different phases, and hence 
a high entropy S. In order for S to remain constant as 
wc move away from the boundary, the temperature must 
increase. If an experiment were performed in which D 
were ramped, the lattice would cool as the phase bound- 
ary is approached from the ferromagnetic side, and heat 
as one moves away into the paramagnet. 

Comparing the isentrope contours, Fig. [4j with the 
those for spin density as shown in Fig. El we se e that as we 
traverse the first order phase boundary, the spin density 
changes more rapidly as the temperature T approaches 
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FIG. 5: (color online). Ground state MFT phase diagram of 
the two dimensional square lattice Hubbard model. Antifer- 
romagnetism is favored at and near half-filling, and extends 
all the way to U = as a consequence of the divergence of the 
p — 1 density of states. Ferromagnetic regions are present at 
stronger coupling. 



FIG. 6: (color online). Isentropic curves of the half- filled 
square lattice Hubbard model. S is independent of U in the 
paramagnetic phase. The isentropic curves are continuous 
and bend upwards upon entering the ordered antiferromag- 
netic region. The entire T — axis is AF. (See text and 
Fig. El) 



zero. In fact, this is consistent with Clausius-Clapeyron 
equation which relates the slope of the phase boundary 
with the change in entropy and spin density by, 

dT) Sf m Sp m . . 

where Sf m and s pm stand for entropy in the ferromagnetic 
and paramagnetic phases respectively. We have verified 
that our results satisfy this condition quantitatively. 

IV. Isentropic Curves of the Hubbard Hamiltonian 

We now turn to the isentropic curves of the Hubbard 
Hamiltonian. The ground state MFT phase diagram is 
shown in Fig. [5j and consists of paramagnetic P regions 
adjacent to empty and fully occupied fillings (p = 0,2). 
In the center, closer to half-filling, magnetic phases arise, 
with antiferromagnetism AF predominating immediately 
adjacent to p = 1 and ferromagnetism F a bit farther 
away, at sufficiently large U . The phase diagram is sym- 
metric about p = 1, as a consequence of the particle- hole 
symmetry of the Hubbard model on a bipartite lattice 
with near-neighbor hopping. 

Figure [5] shows the phase boundary in the T/t — U/t 
plane at half- filling p=l. As discussed earlier, the asso- 
ciated isentropes are parallel to the U/t axis in the para- 
magnetic phase. They then bend to higher T/t as U/t 
increases in the antiferromagnet. This adiabatic heating 
is explained by the same reasoning as for the classical 



BCM case: The phase boundary represents a location of 
particularly high entropy, so that for S to remain fixed 
as one leaves its vicinity the temperature must increase. 
The gross morphology for the isentropes found in Fig. [6] 
agrees well with exact QMC calculations^ which can be 
performed with no sign problem in this half-filled case. 
The QMC method, which is exact, of course captures the 
fact that there is no finite temperature phase transition 
in the square lattice FHM. The role of the phase bound- 
ary there is played by the temperature scale of the mean 
field charge gap, which is determined by the plateau in p 
versus p. 

Figure [7J shows the analogous phase boundaries and 
isentropes for p — 0.8. Here the structure is much richer 
since there are three possible phases for this density, with 
the paramagnet first giving way to antiferromagnetic or- 
der as U /t is increased, followed by a second transition 
to ferromagnetism. In this case the isentropes can bend 
either to higher or to lower T as the AF boundary is left 
with increasing U . The decrease in T occurs at lower T/t 
where the F boundary is more proximate to the AF one. 
The isentropes move to higher T as U/t decreases from 
the F phase boundary. The other interesting feature of 
Fig. [Jj not present in Fig. [6j is the discontinuity in the 
isentropes at the AF-F boundary. This occurs here, sim- 
ilar to the situation for temperatures below the tricritical 
point in the BCM model, because of the first-order na- 
ture of the transition. As with the BCM, we have verified 
that the entropy jump along the first order boundary sat- 
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FIG. 7: (color online). Isentropic curves of the square lattice 
Hubbard model at density p — 0.8. As at p — 1.0 (Fig. |SJ, 
S is independent of U in the paramagnetic phase. The isen- 
tropic curves are continuous upwards upon entering the or- 
dered antiferromagnetic region, trending upwards for large S 
and downwards for small S. The isentropic curves are dis- 
continuous at the AF — ► F transition. As U is increased at 
T = 0, there is first a P — > AF transition, followed by a AF 
—J- F transition. (See text and Fig. [5]) 



isfies the Clausius-Clapeyron relation, which here takes 
the form, 



dU 
dT 



Sfn 



^afm 



fl.fi 11 



(10) 



Here ^ is the slope of the AF-F phase boundary, 
Sfm,Safm are the associated entropies, and 2?f m ,2? a fm are 
the double occupations in the two phases. 

The behavior of the ferromagnetic, mf, and antiferro- 
magnetic, m a f, order parameters for p = 0.8 and T = 0.2 
is shown in Fig. [8] Both change discontinuously through 
the first order AF-F transition at U rj 8.8t. However 
at the second order P-AF transition at U « 3i, m a f in- 
creases continuously from zero with the MFT exponent 

An interesting feature of the P — > AF phase bound- 
ary at p = 0.8 is that the critical U initially decreases as 
the temperature T increases. Put another way, as T is 
lowered for U « 3.2 there is a P — > AF ordering transi- 
tion, but then there is a rc-cntrance to the P phase as 
T is reduced further. We have verified that this phe- 
nomenon is obtained also within an independent random 
phase approximation calculation, in which the magnetic 
susceptibility is given by, 



x(q,T) 



xo(q,r) 

l-C/ X o(q,T) 



X0 (q,T) = 



/(CM 



/(eu) 



£k — e k+q 



(11) 



The critical temperature is then determined from Eq. [TT] 
via the Stoner criterion 1 — t/xo(q, T) = 0, and exhibits 
the same re-entrant phenomenon as the MFT calculation. 

Finally, we show in Fig. [9] the phase diagram and isen- 
tropes for quarter filling. At zero temperature (see Fig. [5]) 
there is a paramagnetic to ferromagnetic transition as 
U increases. However, at higher T an intermediate an- 
tiferromagnetic phase intervenes. The general trend is 
towards adiabatic heating as U rises. 



V. Conclusions 

In this paper, the isentropic curves of a classical sys- 
tem of magnetism, the Blume-Capel model, and of the 
quantum fermion Hubbard Hamiltonian, have been de- 
termined. 

In the case of the BCM, Mean Field Theory and Monte 
Carlo (Wang-Landau) calculations give a qualitatively 
similar pictures in which adiabatic heating is observed 
as one moves away from the phase boundary, although 
the precise, quantitative location of the transition is, of 
course, different in the two methods. The behavior of 
the isentropes is made somewhat more complex by the 
presence of a tricritical point on the BCM phase bound- 
ary, so that there is a region at low T and large vacancy 
fraction where the curves are discontinuous. 

For the FHM, we have presented only MFT results, 
since the sign problem in general prevents Quantum 
Monte Carlo simulations at low temperatures. The topol- 
ogy of the isentropes, like the phase diagram itself, is 
relatively simple at half-filling, consisting of lines paral- 
lel to the U axis in the paramagnetic phase, and then 
trending upwards as the antiferromagnetic boundary is 
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FIG. 9: (color online). Phase diagram and isentropes for 
quarter filling. For T > 0.5£, an antiferromagnetic region 
appears between the paramagnetic region at weak coupling 
and ferromagnetic region at strong coupling. 

crossed. For this half-filled case, QMC is possible, and 
we have compared our MFT results to those calculations. 



The isentropes are considerably more elaborate in the 
doped case, where MFT exhibits competing transitions 
to ferro- and antiferro-magnetic orders. In these cases 
the dependence of the temperature along the isentropes 
is non-monotonic. 

It is well known that MFT predicts some qualitatively 
incorrect features of the FHM phase diagram, including, 
for example, the existence of a finite temperature Need 
transition. Nevertheless, the rough morphology of the 
isentropes within MFT and QMC are similar, with the 
role of the MFT phase boundary played by the temper- 
ature of the charge gap, which is clearly seen to open 
at nonzero temperatures. This provides some assurance 
that the MFT results will remain qualitatively accurate 
in the doped case, even in the absence of exact QMC 
results at low temperatures there. 

Optical lattice experiments for bosons 3 -^— and 
fermion o 42 ' 43 have reached temperatures (T ~ t in the 
bosonic case) that provide strong evidence of the Mott 
transition. A reduction in the experimentally accessible 
entropy per fermionic atom, s ~ ln2, by a factor of 2-3 
will allow the observation of local spin correlations^ One 
possibility raised by the results of Fig. [7] is that adiabatic 
cooling can occur in the proximity of two competing types 
of order. However, Fig. [7] also suggests that this cooling 
only occurs when T is already sufficiently low. 

We acknowledge financial support from ARO Award 
W911NF0710576 with funds from the DARPA OLE pro- 
gram. 
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